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> ' Abstract 

r» : 

Cf") , We investigate the formation cross sections of the deeply bound 

pionic atoms on halo nuclei theoretically. We consider one-neutron 

■ pick-up (d, 3 He) reaction for the pionic atom formation and find that 

■ the formation rate is significantly increased by the existence of the 
halo neutron with long tail of the wave function. This enhancement 
is expected to increase the experimental feasibility for the formation 
of the deeply bound pionic states on /3-unstable nuclei. 

'O 

1 Introduction 

> ■ 

^ ■ Mesic atoms are the useful laboratories for studying the meson-baryon inter- 

actions, the meson properties in nuclear medium and the nuclear properties. 
Among others, pionic atoms have been long utilized for such studies [1]. The 
binding energies and widths of the pionic bound states are expected to pro- 
vide us with the unique information of the pion-nuclear interaction [2] and 
of the nuclear density distributions [3, 4]. 

The standard method to produce pionic atoms starts with injecting slow 
negative pions into matter. The pions will be stopped and trapped in out- 
ermost orbits of atoms, then lose energy by emitting Auger electrons and 
x-rays to cascade down to deeper atomic states. However, the x-ray cascade 
ceases at the so called 'last orbital', where the pions are absorbed by the nu- 
cleus without going to deeper atomic states. This prevents us from studying 
the deep pionic states for heavy nuclei that have large overlap with nuclear 
density [5]. 

To overcome this difficulty, the recoilless mesic atom production by (d, 3 He) 
reaction has been proposed [6]. This method tries to directly create deeply 
bound pions by the nuclear reaction without cascade down of the pion from 
the outer most orbits. This direct method was proved to be effective and 
powerful for the study of far deeply bound pionic atoms [7, 8]. We can now 
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determine the binding energies and widths of the deepest bound states such 
as Is and 2p atomic states from the energy spectra of the emitted 3 He. We 
remark that the (d, 3 He) reaction is so powerful that it may also be applied 
to the creation of neutral meson-nucleus bound states such as the i] and u 
mesic nuclei [9, 10]. 

On the other hand, the physics of the unstable nuclei is one of the most 
interesting subject in modern nuclear science. Since the successful use of 
a new experimental tool, namely "beams of unstable nuclei" [11, 12], the 
properties of nuclei far from the stability line have been studied extensively 
in many laboratories where secondary beams of unstable nuclei are available. 
One of the exciting findings was the neutron halo structure around the 9 Li 
core in n Li [11]. Since then, the structure and reactions of /^-unstable nuclei 
have been studied extensively both theoretically and experimentally [13, 14]. 

Since the repulsive S-wave pion-nucleus interaction causes the formation 
of the pionic halo around the core nucleus [5], the binding energies and widths 
of deeply bound pionic atoms are expected to be very sensitive to the struc- 
ture of the nuclear surface, in particular, to the neutron skin [4]. Hence, 
we can expect to obtain new information on neutron/proton distribution of 
the unstable nuclei by observing the deeply bound pionic atoms. The deter- 
mination of the neutron distributions for unstable nuclei is one of the most 
important subjects in the /5-unstable nuclear physics [13]. In addition the 
(d, 3 He) reaction with pionic atom formation proceeds through the single- 
neutron pick up and are known to be very sensitive to the single-particle 
properties of the neutron in nucleus [15]. Hence, we can also expect to 
obtain the new information of the neutron single-particle states (nlj) in un- 
stable nuclei. This is very interesting in the context of the possible change 
of the magic numbers far from stability line [16]. 

For the formation of the deeply bound pionic atoms on unstable nuclei, 
we use the inverse kinematics technique, in which a light target such as 
deuteron is bombarded by radioactive ion (RI) beam of intermediate energy 
and the recoiling light ejectile 3 He is detected after forming pionic atoms 
in the heavy projectile [17]. We evaluate the formation cross sections for 
inverse kinematics reactions d(RI, 3 He)X in order to examine the experimental 
feasibility [18]. Recently, the existence of the halo like structure of neutron 
distribution is predicted theoretically in medium heavy mass nuclei like Zr 
isotopes [19], which was not considered in our previous work [18]. This is very 
interesting since the existence of the long tail of the neutron wave function 
may enhance the formation cross section of the deeply bound pionic atoms 
in unstable nuclei and increase much the experimental feasibility in very 
unstable region of nuclear chart using the next generation of accelerators 
such as RIKEN RIB factory. 

In this paper we study the formation cross sections of the deeply bound 
pionic states for the medium and heavy halo nuclei. We simulate the halo-like 
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structure by using the solutions of the Schrodinger equation with an one-body 
neutron potential for various potential depths, and we investigate the effects 
of the long-tail of the neutron wave functions to the (d, 3 He) spectra. In 
section 2, we describe our theoretical formalism. We show the numerical 
results in section 3. Section 4 is devoted to the summary. 



2 Formalism 

In this section, we describe our theoretical formalism to investigate the pionic 
atom formation spectra on halo nuclei. As we will describe below, we intro- 
duce a single free parameter F, a reduction factor for neutron potential depth, 
and simulate the halo structure by changing the parameter F. Nuclear prop- 
erties, such as neutron wave functions, density distribution and one-neutron 
separation energies, are obtained consistently for various F values. Pionic 
atom structure and formation spectra are also calculated consistently for 
several F values. The factor F is the only free parameter introduced in the 
present work. 

In order to simulate the halo wave functions for neutron, we solve the 
one-body Schrodinger equation for all occupied single-particle states for both 
protons and neutrons using the potential in ref. [20]. The one-body potential 
for protons and neutrons are given by, 



1 + exp ((r - Rq) /a ) 
h 2 . . 1 d ( V \ wi\ 

~2MV ( ' Sj r ~dr \ 1 + exp ((r - R ls ) /a ls ) ) + W ' ( j 

where V c is the Coulomb potential of a uniform charge distribution with the 
radius R c . We adopt the potential parameters which were determined so as 
to fit the experimental single particle energies for Pb region in ref. [20] and 
these are listed in Table 1. 

Here, we introduce a free parameter F which is a reduction factor for 
neutron potential depth to simulate the halo structure as; 

V - F x V . (2) 

This factor is introduced only for the neutron potential and the potential 
depth for proton is not altered for all numerical calculations in the present 
work. 

Using the proton and neutron radial wave functions R^ij^ ( r ) , we calculate 
the proton and neutron density distributions by summing up the absolute 
square of the wave functions for all occupied states (nlj) for proton and 
neutron, respectively, 
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Pip/n) (r) 
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(3) 



We calculate the density of the target nucleus pa using Eq.(3), and we eval- 
uate the final-state nuclear density Pa~i, which appears in Eq.(14), by re- 
moving one neutron in the state which has the largest contribution to the 
(d, 3 He) spectra from the pa- 

The pionic atom wave functions and eigen energies are obtained by solving 
the Klein-Gordon equation with pion-nucleus optical potential, 

-V 2 + p 2 + 2pV opt (r)] ^ (r) — [E — V Co ui (r)f 0. (r), (4) 

where p is the reduced mass of pion-nucleus system and Vcouiij') is the 
Coulomb potential with a finite nuclear charge distribution p p (r) given by, 



Vcoui (r) 
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Ir — r'l 



d 3 r'. 



(5) 



We use the proton and neutron density distributions p p and p n of target nu- 
cleus given in Eq. (3) for all densities in the Klein-Gordon equation. Hence, 
all nuclear densities and the pion wave function depend on the parameter F. 

For the pion-nucleus optical potential, we adopt the conventional Ericson- 
Ericson form [21]; 



2pV opt (r) = -4tt b (r) + e 2 B p 2 (r) 



with 



4ttV • 



c(r)+ejV(r) L(r)V, 



(6) 



b (r) = ei [b p (r) + b x [p n (r) - p p (r)]], 
c (r) = q 1 [c p (r) + ci [p n (r) - p p (r)]], 



i(r) 



4 x 
1 + -.A 



c(r) + e 2 L C p z (r) 



(7) 
(8) 

(9) 



We adopt a parameter set determined by Seki and Masutani in ref. [22] as 
the standard one and these are listed in Table 2. 

To evaluate the formation rate of the pionic atoms in the (d, 3 He) reac- 
tions, we use the effective number approach [15] which is known to predict 
the experimental cross sections reasonably well [7]. The (d, 3 He) spectrum is 
given as; 



d 2 a 



dQuedEne , 



dA- > 3 He(A-l) 7 r 



i,W^ cff 2vrA^ + rL/4' 



(10) 
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with 



AE = Q + 77v - BE* + S n - [M d + M n - M He ], 



(11) 



where BE n and are the binding energy and the width of pionic atom, 
and S n is the neutron separation energy. These quantities also depend on 
the factor F introduced in Eq.(2). 

The effective number N c g in Eq.(lO) is defined with the neutron and the 
pion wave functions ip n (r) and <f> n (r) as, 



N eS = J2 I X) (r) [C (0 ® ^ (r)] JM Xl (r) d\ 

JM J 



(12) 



where the deuteron and 3 He wave functions are expressed as Xi an d Xfi anc l 
are calculated in Eikonal approximation as; 

X/( r )^(r) =exp(iq-r)D(z,h). (13) 
The distortion factor D (z, b) appeared in Eq. (13) is defined as; 



D (z, b) = exp 



1 r z 1 f°° 

-^ a dN / dz' p A (z, b) - -a hN / dz' p A -i (z', b) 



(14) 

where the nuclear densities are given in Eq. (3). 

We like to mention here again that we have only a single free parameter F 
in this model and all physical quantities are obtained in a consistent manner 
for each F value to simulate the effects of the halo structure on the pionic 
atom formation (d, 3 He) spectra. 





Vo(MeV) 


i?o(fm) 


cco(fm) 


Ms 


i^(fm) 


a is (fm) R c (im) 


protons 


-60.94 


7.52 


0.79 


22.33 


7.28 


0.59 6.70 


neutrons 


-49.40 


7.05 


0.66 


21.47 


6.78 


0.24 



Table 1: Parameters for proton and neutron one-body potential obtained in 
ref. [20] for Pb region. 



3 Numerical results 

In this section we show the numerical results for the pionic atom formation 
spectra in (d, 3 He) reaction by varying the parameter F to simulate the halo- 
like structure. First, we consider the pionic atoms on 208 Pb. The calculated 
neutron binding energies are listed in Table 3. Obviously, the binding energies 
are smaller for smaller F values. The H13/2 state does not bound for the 
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b = - 0.0283m; 1 61= - 0.12m; 1 
c =0.223m; 3 ci=0.25m; 3 
^0=0.042 i m~ A C =0.10 i m~ 6 
A=1.0 

Table 2: Pion-nucleus potential parameters obtained in ref. [22]. 



values F smaller than 0.765. Thus, we vary the F between 1 and 0.765. These 
binding energies are used as the neutron separation energies S n appeared in 
Eq.(H). 

We show the nuclear density distributions in Fig. 1 for the values F — 1.0 
and 0.765. We can see that the neutron density distribution is affected signif- 
icantly by the choice of F value and has longer tail for F=0.765. The nuclear 
densities shown in Fig. 1 are used to calculate the pionic wave functions. 

The density distribution of the neutron 3pi/ 2 state and that of the pionic 
Is state are shown in Fig. 2. We can see that the neutron single-particle 
density in the valence level is sensitive to the F value and has halo-like tail 
for the smaller F values. On the other hand, the density of the pionic atom 
depends on the F value in rather mild manner as seen in Fig. 2. The binding 
energy of the Is pionic state is 6.93MeV for F—l while 6.67MeV for F=0.765, 
respectively. 

We also show the calculated root mean square radii for valence neutron 
levels and show the dependence on the F value in Table 4. We can see that the 
neutron states with lower angular momenta are more sensitive to the F value 
and develop the longer tail. This feature is interpreted as the consequence 
of the centrifugal barrier and was suggested in ref. [23] . We expect that the 
long tail of the neutron distribution leads to larger cross section for producing 
deeply bound pionic states since the quasi-substitutional states are know to 
be largely populated for the recoilless (d, 3 He) reactions [6]. 

To see this, we show the calculated (d, 3 He) spectra for 208 Pb target in 
Fig. 3 as a function of Q-value. Here, we only evaluate the contributions 
from bound pionic states. As can be seen in the figure, we find that the for- 
mation cross section is significantly enhanced for the case of halo-like neutron 
distribution corresponding to F=0.765. 

The formation cross section of [(2p 3 / 2 ); 1 (X>(2p) 7r ] configuration, which has 
the largest contribution to the spectra, is around factor 3.5 larger for F=0.765 
than that for F=1.0 which corresponds to no-halo case. We also find that 
the shape of the spectrum is considerably influenced with the change of the 
F value. This is due to (1) the change of the separation energy of each 
neutron state as shown in the Table 3, and (2) the change of the momen- 
tum transfer for different Q-values. Thus, for example, we can see from 
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Fig. 3 that the energy difference between two peaks [(2p3/2)^ 1 ®(2p) 7r ] and 
[(2p 1 / 2 ) _1( S>(2p) 7r ] are much closer for F=0.765 case and there appears the 
extra strength around Q=-127 ~ -128 MeV due to the [(2p 3/2 ) T 7 1 ®(ls) 7r ] and 
[(2p 1 / 2 )^ 1 ®(ls) 7r ] contributions, which are negligible for F=1.0 case because 
of the smaller momentum transfer. 

We also investigate the case of Xe, which is known to have large formation 
cross section for pionic Is state coupled to the Si/ 2 neutron state [15]. For 
this case, we used the radius parameters, i?o, Ris and R c , for the proton and 
neutron potential listed in Table 1 by assuming the mass number dependence 
A 1 / 3 . The other parameters are assumed to be the same. We show the 
calculated binding energies in Table 5. We do not have the bound lhn/ 2 
state for the value of F smaller than 0.775. The calculated root mean square 
radii are shown in Table 6. We find again that the long halo-like tail is 
developed well for 3si/ 2 state, the state with the lowest angular momentum, 
as in the case of Pb. 

We show the calculated spectra for the case of Xe target in Fig. 4 as a 
function of the reaction Q-value. We can see again that there exists large 
enhancement of the cross section for halo-like neutron states. We also find 
the change of the spectrum as in the case of Pb. The formation cross sections 
of [(3si/ 2 )~ 1 ®(ls) 7r ] configuration, which has the largest contribution to the 
spectra, is around factor 2.9 larger for F=0.775 than that for F=1.0 which 
corresponds to no-halo case. 

Finally, we examined the sensitivity of our results to the theoretical in- 
puts. In Fig. 5, we show the (d, 3 He) spectra for Pb which are calculated by 
using the nuclear density with the relativistic mean field model (RMF) [24] 
and with the pion- nucleus optical potential obtained by Konijn et al. [25]. 
The RMF model gives the nuclear density without the halo-like structure. 
For comparison, we also show the results calculated by the nuclear density 
with the potential parameters in Table 1 denoted as p spe th [20] and/or the re- 
sults with the pion-nucleus optical potential by Seki and Masutani [22] (SM). 
As seen, the calculated spectra are slightly modified by the above choice of 
the nuclear densities and the pion-nucleus optical potential. But these modi- 
fications are not significant for the cases considered here. Next, in Fig. 6, we 
show the results by assuming the halo-like nuclear densities. We adopted the 
same procedure as in Section 2 to simulate the possible halo-like structure 
of the nuclear density. The strength of the single-nucleon potential is varied 
by introducing a single free parameter F as in Eq.(2). The nucleon single- 
particle potential by Koura (set OB) [26] is adopted to see the dependence 
on the nucleon-nucleus potential. For this potential, we found F=0.715 is 
the smallest F value which gives the bound lii3/ 2 single-particle state. Since 
the overlap of the pion wave function and the nuclear density is larger for 
this case, we obtained significantly larger widths for the deeply bound pionic 
states. Thus, the calculated spectrum with the nuclear density by the single 
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nucleoli potential by Koura is influenced as seen in Fig. 6: (i) the largest 
two peaks due to [(2p 3 / 2 )^ 1 ®(2p) 7r ] and [(2p 1 / 2 )^ 1 ®(2p) 7r ] configurations can 
not be distinguished, and (ii) the bump at Q=-127 ~ -128MeV due to pionic 
ls-states contribution is not seen with the nuclear density by Koura. How- 
ever, we can see that the absolute heights of the peaks in the spectra are 
roughly the same. We can conclude that, though the detailed shape of the 
spectrum is influenced by the theoretical inputs, the formation cross section 
of the deeply bound pionic states is largely enhanced by the existence of the 
halo-like structure in the target nucleus. 



B.E.(keV) 


3pi/ 2 


2f 5/2 


3p 3 /2 


H13/2 


2f 7/2 


lh 9/2 


F=1.0 


7761.6 


8386.0 


8807.6 


8895.0 


11252.3 


11868.5 


F=0.765 


1094.4 


933.91 


1513.8 


18.281 


2787.6 


2864.4 



Table 3: Calculated neutron binding energies for each valence neutron level 
of Pb in unit of keV for F=1.0 and 0.765. 



< r 2 > i /^(f m ) 


3pi/ 2 


2f5/2 


3P3/2 


lil3/2 


2f7/2 


lh 9 /2 


F=1.0 


5.9145 


5.7722 


5.7723 


6.2341 


5.6795 


5.6008 


F=0.765 


7.9384 


6.9139 


7.4099 


6.5740 


6.3475 


5.9810 



Table 4: Calculated root mean square radii for each valence neutron state of 
Pb in unit of fm for F=1.0 and 0.765. 



B.E.(keV) 


lhn/2 


3si /2 


2d 3/2 


lg7/2 


2d 5/2 


F=1.0 


8215.9 


9668.8 


9219.1 


11836.0 


11726.4 


F=0.775 


121.0 


2821.6 


2306.4 


3579.1 


3946.9 



Table 5: Calculated neutron binding energies for each valence neutron level 
of Xe in unit of keV for F=1.0 and 0.775. 



4 Conclusions 

In this paper, we have studied theoretically the formation cross sections of the 
deeply bound pionic states in (d, 3 He) reactions for halo nuclei. To simulate 
the halo-like structure of neutrons, we have simply modified the depth of the 
neutron potential by introducing a single free parameter F and generated the 
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< r 2 > 1/2 (fm) 


lhn/2 


3si/a 


2d 3/2 


lg7/2 


2d 5/2 


F=1.0 


5.4817 


5.1079 


5.0575 


4.8649 


4.9911 


F=0.775 


5.8531 


6.2414 


6.0258 


5.2500 


5.614 



Table 6: Calculated root mean square radii for each valence neutron state of 
Xe in unit of fm for F=1.0 and 0.775. 



0.12 



0.08 



0.06 




0.04 



0.02 



Figure 1: Neutron density distributions p n are plotted as a function of the 
radial coordinate r for F=1.0 and 0.765. Proton density p p is also plotted as 
the solid line. 



shallow neutron states with long tail. Using the neutron wave function and 
the nuclear densities thus obtained, we have calculated the formation cross 
section of the deeply bound pionic states in (d, 3 He) reactions. By varying the 
parameter F, we can see the effect of valence neutron states on the formation 
spectra. We found that the formation cross section is significantly enhanced 
for the case with neutron halo-like states. We have examined the several 
theoretical models and found that the enhancement of the reaction strength 
due to the halo-like structure in the target nucleus is the common feature 
of the spectrum. For the cases with halo-like states, the reaction strengths 
are enhanced about a factor 3.5 for Pb and about 2.9 for Xe than the case 
without neutron halo. We also found that the shape of the reaction spectra 
is considerably affected by the choice of the strength parameter F for the 
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Figure 2: Density distributions of the neutron 3p!/ 2 and the pionic Is states 
are shown for F—l.Q and 0.765 for Pb. 



nucleon-nucleus potential. 

We believe that the present results are relevant to the development of 
spectroscopic studies for the deeply bound pionic states in wide area in nu- 
clear chart including the neutron-rich unstable nuclei. 
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